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Outline 


- Ionization streamer formation and propagation synopsis 

- “NSLE” propagation ionization numerical simulations 

- Reported data and theory 

- Analysis of the possibility of radiated fields 
and corroborating Model 

- Maxwell numerical simulations 
and test case results 

- Conclusion and comments 








RF Radiation from Laser Filaments 


• THz and sub-THz radiation from short pulse laser filaments have been 
reported. 

• Sub-THz radiation from longitudinal plasma filament oscillations has 
been predicted observations reported by Tzortzakis et al (02). 

- Cheng et al (01) have proposed plasma oscillation mechanism 

- This mechanism has been challenged (Carron, 02) 

— Carron states filament can only radiate from discontinuities when filament 
starts and stops. 

- Traveling wave model supports Canon’s results (Zimmerman, 02) 

• At present there is a shortage of quantitative data and no satisfactory 
theory for observed sub-THz radiation 


Comments on this Chart: 

Among the interesting phenomena generated by these laser filaments in air 
is the emission of THz and sub THz electromagnetic radiation. Both near 
and far field measurements have been made by various investigators. 
Unfortunately the reported data is qualitative in nature. A model of sub THz 
radiation was proposed by Cheng et al in which longitudinal oscillations are 
excited by the laser pulse. This model has been challenged on two grounds 
The first is that a traveling current pulse is not expected to radiate except at 
discontinuities. The second is that the electron neutral collision frequency is 
on the order of 10 A 12 sec-1 and so sub THz oscillations would be expected 
to be damped out. At present there does not appear to be a satisfactory 
model to explain the experimental results 








Analysis of Filament Fields 



Phase 1: Initial Dipole !s Created 
•“Capacitor" Is Created 

•Electrons “Pushed" Ahead by 
•Radiation Pressure 
•Ponderomotlve Force 

•Changing Dipole Moment Creates Radiated Field 



Phase 2: Dipole Moves with Laser Pul se Leaving Plasma Trial 
•Two Possibilities 

♦“Capacitor* Moves with Uniform Velocity & 

Trailing Neutral Plasma 
•Electrons Move with Uniform Velocity & 

Trailing Positive Plasma 
•Neither Case Creates Radiation 



Phase 3: Dipole Decays 

■Both Phase 2 Possibilities Finally Leave Neutral Plasma 
•Dipole Creates Radiation as it Dim inishes 


Comments on this Chart: 

1 . initial dipole is created, a radiated field is created since there is a 2 nd time 
derivative of the dipole moment 

2. The pulse propagates, there are only linear rates of dipole moment change 
so the 2 nd derivative of P is zero and there is no radiated field 

3. The dipole is extinguished, a radiated field is produced since there is a 2 nd 
time derivative of the dipole moment 






Diagram of Possible Ionization Streamer Oscillations 
from Scaled Plasma Frequency Parameters 



Wp is the plasma 
frequency 

V is the collision 
frequency 

a is the channel 
radius • 








Comments on this Chart: 

This shows the derivation path from Maxwell’s Equations to the Envelope 
Equation 

1. With Uniform Media and no free charge, Maxwell’s equations becomes the 
wave equation 

2. Envelope approximation, expansion of linear polarization in terms of time 
derivatives of the envelope and transformation to retarded time gives final 
envelope equation 







Real Time Maxwell Algorithm 


Instead of the using a wave equ ation 
the curl equations are directly solved 
simultaneously 

p 

VxE=-B 

□ 

Vx/f = -D+J/ 


The media polarization current is 
determined from a nonlinear 
oscillator model analyzed (via ODE’s) 
in conjunction with the field time 
stepping 


Using the Yee algorithm the E field and H field 
meshes are offset in space allowing centered 
2nd order curl differencing for finding E or B 
at a new time step from H or E at the o Id time 



Comments on this Chart:. 

Start with Maxwell’s equations 

Yee Algorithm has Electric Field and Magnetic Field offset 
in space and time 











Inclusion of Frequency Dependant Properties 
In The Time Domain Calculation 


Non nonlinear media are 
included by solving a simple 
oscillator equation including Linear 
both linear and nonlinear Res P Qn! 
terms 

Simultaneously carrying this 
solution provides a 
polarization current for 
Maxwell’s Eqs. 

The oscillator coefficients are 
determined from optical data 


Displncemen 

t I 


Maxwell Equation 

Vx//=^aD/af+y 7r 


Non Linear (Kerr effect) 


Oscillator equation of motion 

+—(ff* vp-H r oe)j 


Similar terms Tory 
and 2 


-HH^***) 


Polarization Current 


Similar tenns forx 
andy 


Comments on this Chart: 


Solution of Fields are combined with Source term computation 
Shows details of mapping between a harmonic oscillator model and th 
nonlinear effects 

Oscillator model 

1. Begins with linear part 

2. Extended to non-linear (Kerr) effect 

3. Finally results in ionized atom 
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Finite Difference Model ODE’s 


Electric field ODEs for the elemental cell (include corresp onding Y and Z Terms) 

\ E I J \ £ i J 

Dielectric/polarization model, x, — (v y ). (Include corresponding Yand Z Terms) 

(v‘), 4jrX E >)Aj?’i\{ < H - >< M, - <H ^ ('’-) 1 )-(2r)(v,) 1 -(•?)*+«? 


Dielectric/polarization model, n i = S t (Include corresponding Y and Z Terms) 

00,■ = >/ OO,"^ 00,) + [^)00,”, 


Ionization model 


S(e/cm>-s e c)=^-=»(W 2 )* w , (///,r 


Comments on this Chart: 

This shows how the physics equations become algebraic equations to be 
solved 

1 The electric field depends on the linear and non-linear sources 

2. Bound electrons move under the effects of the electric and magnetic fields 
as well as the restoring forces of the oscillator model for bound electrons 

3 f ree electrons move under the effects of the electric and magnetic fields. 

4. The equation for the ionization is also shown 





Computational Check out problem 


A femtosecond pulse propagating in air impinges 
on a Silica slab. 


This problem has been calculated to check out the 
Maxwell formulation and results. 



The form of the incident pulse in the 
these code checkout calculations is 
taken to be: 

E(t) = E„ exp f ~ ~ jsin(2ff/ 0 (f-r)) 


T = 8.333/5, / 0 = 0.3 pHz , r = 33.33/s 








Field Propagation Checkout -1 















Field Propagation Checkout Results - II 


Willi Cubic Term . W/O Cubic Tcnn — 


1“ Kerr effect on speed in Si (20 GV/M) 

2- Induced current in air plasma at lower intensity 

3" Induced current In air plasma at higher intensify 
Showing non linear behavior 



i 
























Comments on this Chart: 


The question of whether RF radiation can be generated from a traveling 

Current pulse can be found by integratin over the moments of the dipole 
distribution, accounting for time retardation (coherence). 


The usual expansion of the fields into near intermediate and radiated 
components can be found. 






An alternate derivation is to proceed from the propagating current distribution. 

Standard expressions for the radiated field can be evaluated. This of course gives the 
same result as found from the dipole moment integration. 
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Conclusions 


I - First principles analysis indicates no radiated field from a uniformly 
propagating ionization streamer 

2- The predicted ionization densities/collision frequencies would 
abruptly damp out any oscillations in the ionization 

3- Maxwell simulations showed both optical frequencies and a low frequency 
component. The low frequency is not low enough to explain reported sub THz 
radiation. 

4- Observation show numerous streamers in one laser beam. 

NLSE simulations indicate the formation and extinction of multiple ionization 
Streamers. These non steady state features have not been analyzed here 

5- Our simulation do not indicate any sub THz radiation. 

6- Mechanisms to explain reported observations are not well 
understood. 







Test Case, Dielectric Model (Continued) 


Analytic Estimate for the transverse and longitudinal Current Densities 


transverse 

z(polarization) = 7.02x10" 12 tw 
v z (polization) ~ 1.32 x 10 4 m ts 
y,s5,61xlO w >l/« a 


longitudinal 

x(polarization) si.55x10" 56 m 
v T ( polarization ) = 5.48 x 10" 1 m / s 
J x ~2A^x\0 6 Alm 2 


Numerical results for the transverse and longitudinal current densities 

transverse . longitudinal 

z{polarization)^\.2\y.\(y n m x(polarization) ~ 4.65xl0" l6 m 

v .{polarization) s 2.25xl0 4 /n/s v x (polarization) s 1.55 x10 I) /w/j 

J s = 9.40x1 O lo A/m 2 J x s7,2Axl0 6 A/m 2 

The analytic and numerical results are in reasonable agreement 


Analytical and numerical results are shown on this slide. 


As can be seen the analytical and numerical results are in good agreement, 
comparison with the previous slide 
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